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The effect of limited depth 
and width of waterway on performance 
of ducted propellers 


Wojciech Gorski 
Ship Design and Research Centre S.A., Gdansk 


Jan Kulczyk, Tomasz Tabaczek 
Wroctaw University of Technology 


ABSTRACT 


Model tests of propeller performance in bollard conditions, in deep and shallow water, were carried out at 
Ship Design and Research Centre in Gdansk. Corresponding calculations of propeller performance with ac- 
count for finite dimensions of canal cross-section were carried out at Wrocław University of Technology by 
using their own theoretical model of propeller — hull interaction. The calculations were carried out in model 
scale, at the same water depth as in model tests. For given hull form, propeller geometry and canal cross- 
-section the HPSDK computer code was used to calculate wake fraction, as well as propeller thrust, torque 
and efficiency. The distribution of pressure on waterway bottom and ship sinkage were also determined. 


Keywords : inland navigation, performance of ducted propeller 


INTRODUCTION 


The effect of limited depth and width of waterway on the 
performance of propellers, as is the case in inland navigation, 
have not been by now extensively reported in the literature. In 
the RTD project INBAT the above problem was considered 
important, for it is related to efficiency of propulsion. In order 
to investigate mutual interaction of ship and waterway both 
experimental and theoretical research tools were involved. 
Because of confinements imposed by actual diameters of stock 
propellers used in tests and large dimensions of model ship the 
model tests in shallow water were limited to bollard conditions 
(zero speed). Theoretical CFD methods do not impose limits 
on dimensions. However, many simplifications applied to the 
computational method, both to geometry and governing equa- 
tions, make the applicability of those methods restricted. In the 
present investigations the experimental and theoretical methods 
were applied to the same ship at the same operating condi- 
tions, and provided with complementary information. Theore- 
tical calculations were carried out in model scale in order to 
avoid extrapolation errors when analyzing results. 


HPSDK COMPUTER CODE 


When developing the HPSDK code it was assumed that 
hull, propeller (either screw propeller or ducted propeller) and 
waterway bed (banks and bottom) are considered separately, 
and their interaction is taken into account by using iterative 
solution of flow including updated disturbances. That appro- 
ach was forced by low computational power of available hard- 
ware. 3-D potential flow bounded by ship hull, rigid water sur- 
face and canal bed is calculated using the technique of surface 


vorticity distribution. Viscous effect on velocity distribution in 
propeller disk is taken into account by simplified calculation 
of viscous flow between hull and canal bottom. The Vortex 
Lattice Method was applied to propeller flow at a given distri- 
bution of inflow velocity. The above mentioned methods are 
described in details in [1] and [2]. The flow chart of calcula- 
tions is presented in Fig. 1. 


The computations are organized in 5 steps : 


Step 1 

Within this step the necessary input data are prepared for 
calculations. Hull form is defined and processed in SIATKA 
module. The data file is next read by the HPSEDKDB data 
editor that allows to enter a number of complementary data 
such as : water depth, dimensions of canal cross-section, ship 
speed, propeller diameter and the number of vortex panels re- 
presenting water surface around the ship. The HPSEDSRC 
editor is dedicated to propeller geometry, and the HPSEDDSA 
editor to nozzle geometry in the case of ducted propeller.On 
the basis of the entered hull form the VORNET module com- 
putes the coordinates of grid representing the hull in subsequ- 
ent calculations. The TARCIE module calculates the effect of 
viscosity on distribution of flow velocity under the hull. 


Step 2 


The HDKDK and HDKRK modules calculate potential flow 
around the ship and wake fraction, respectively. If the goal of 
calculations is to determine the nominal flow (excluding the 
effect of running propeller), the HDKRK module computes 
pressure distribution on the canal bottom, then ship sinkage is 
determined and the calculations are terminated. If the goal is to 
determine the effective flow and performance of propeller the 
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DATA INPUT 
- Geometry of propeller 


- Hull form 
- Cross-section of canal 


TARCIE 


Potential flow around the hull 


Velocity in propeller disk 
HDKRK 


Flow without propeller 


- Sinkage 
- Trim 


Vortex elements 


representing screw propeller 


Disturbances induced by hull 
in propeller location 


Influence coefficients 
of screw propeller 


HPSGK 


Flow around screw propeller 
HPSWK 


Disturbances induced by propeller 


on the hull HPKDK 


Velocities in propeller disk 
and in other given location 
HPDEF 


Fig. 1. Flow chart of computations carried out by the HPSDK code 
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Y - Flow w/o propeller 


K 


- > 
Gy ~ 


W(K)=W(K-1 


RESULTS: 
- Nominal wake 
- Effective wake 
- Sinkage 
- Trim 


iterative procedure must be completed and the remaining steps 
are pursued. The calculated flow field including viscous ef- 
fects is used as input to propeller flow calculation in the next 
steps. 


Step 3 


Propeller data are prepared for solving the flow around pro- 
peller (the modules : GEOMK, HSMDK, HPDSK, HPSGK, 
HPDGK). A grid of vortex elements representing propeller 
blades and nozzle is defined and the coefficients of governing 
equations are computed. The effects of hull and waterway bot- 
tom are taken into account by calculating the velocity induced 
by vortex elements representing hull and bottom in collocation 
points on propeller. 


Step 4 


Flow around screw propeller is solved by the HPSWK 
module. In the case of a ducted propeller the flow is solved by 
the HPSWK and HPDWK modules running alternately. The 
HPSWK module computes the flow around impeller including 
disturbances (velocities) induced by nozzle. The HPDWK 
module solves the flow around nozzle including disturbances 
induced by the impeller. In the case of a hull-integrated nozzle 
the sector located inside the hull outline is disregarded in com- 
putations. No special adjustment is made in the area of the junc- 
tion. The process of alternate solving the flows around impel- 
ler and nozzle is not illustrated in the flow chart shown in Fig. 1. 
Thrust and torque are calculated by integrating the pressure on 
impeller blades and nozzle. 


Step 5 


Velocities induced by propeller are computed in colloca- 
tion points on hull (by HPKDK module) and in propeller disk 
(by HPDEF module). Effective wake fraction is calculated by 
HDKRK module. 

At this point the convergence of the iterative procedure is 
tested by using the values of effective wake fraction computed 
in the present and previous runs of HDKRK module. If the 
relative difference in wake fraction is less than 0.05 then com- 
putations are terminated. The above convergence value is usu- 


ally reached after 2 or 3 calculation loops including steps 3 
through 5. The values of thrust and torque are averaged becau- 
se the results oscillate. 

Propeller performance, distribution of hull pressure, distri- 
bution of bottom pressure, ship sinkage and velocity distribu- 
tion in flow with and without the effect of running propeller 
are presented below. 


PROPELLER PERFORMANCE 
IN BOLLARD CONDITIONS 


Model tests of bollard pull were carried out at Ship Design 
and Research Centre by using the pushboat model in the scale 
of 1:4.72. The full-scale dimensions of the pushboat were as 
follows : 


Length overall : Loap = 20.5 m 
Moulded beam : Bp= 9.0m 
Moulded draught: dp = 0.6m 


The hull form of the pushboat is shown in Fig.2. The push- 
boat was propelled by three propellers because of the limited 
draught and variable thrust demand depending on a number of 
barges, draught of barges and water depth. Two side propellers 
were the ducted ones with hull-integrated nozzle. The diameter 
of the model propellers used in tests corresponded to 0.690 m 
in full scale. The stock propellers were applied of Ka 4 -70 
Wageningen propeller series, running in 19A nozzle. The cen- 
tral propeller was designed as a ducted propeller of Z-drive 
azimuthing arrangement. It is lowered to operating position at 
a water depth greater than 1.2 m, and retracted when water 
depth is smaller than that. The diameter of the central propeller 
was equal to 1.100 m. In the model tests the CP324 stock pro- 
peller in nozzle 19A was applied [3]. 


Corresponding computations were carried out for the train 
composed of two barges and the pushboat in model scale. Main 
particulars of a single barge were as follows : 


Loa, = 48.75 m 
Bp = 9.00 m 
dp = 140m 


Loap= 20.50 m 
Lwi)= 20.12 m 
Bp = 9.00m 
dp = 0.60m 
V, =82.40 m3 


Fig. 2. Hull form of the pushboat used for model tests and computations 
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Propeller performance 


Results of the model tests and computations are compared 
in Tab.1. The HPSDK code is not capable of computing the 
propeller performance exact in bollard conditions i.e. at zero 
ship speed. In order to overcome the problem the computa- 
tions were run at low ship speed values, namely at 0.75 m/s, 
0.50 m/s, 0.25 m/s, and a constant rotational speed of propel- 
lers. Propeller performance was then extrapolated to zero ship 
speed. The relationship between ship speed and propeller thrust 
appeared almost linear for both central and side propellers (Fig. 
3). The results are consistent with hydrodynamic characteris- 
tics of the propellers. 


Tab. 1. Thrust of the propeller and pull of the pushboat in bollard conditions 


Central 


| 2m | m | 


2S? 75 O 7S || 30 


N 
= 
we 
i=) 
N 
© 


164.9! 72.6 | 173.5! 74.6 


Pull [N]} 218.11 73.8 [208.8; 82.1 [273.5 :109.9 | 296.1; 132.8 


Thrust i 
of i 
impeller ' ' 


IN] 
Tin 211! 76 | 208! 74 | 256! 115 | 260! 116 


0 0.25 


V, [m/s] 0.5 


—@n=12.5rps ,h=2m 
—*h—-n= 7.5rps ,h=2m 


Hn =12.5rps , h=5m 
n= 7.5rps ,h=5m 


Fig. 3. Thrust of the central propeller, extrapolated to bollard conditions 


The following magnitudes were measured 
during the model tests : 


> the net thrust of the entire central propulsion unit 
(including shaft post in front of the propeller) 

> the thrust of the side impellers (excluding the thrust 
of nozzle), and 

> the pull of the pushboat [4]. 


In the computational model the shaft post in front of the 
propeller was not accounted for. The computed thrust of the 
central propeller, shown in Tab.1, contains only the thrust ge- 
nerated by impeller and nozzle. 

Both the results of the model tests and computations reve- 
aled that the effect of water depth on thrust and pull at zero 
ship speed is insignificant. 
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Bottom pressure 


The computations revealed that at the water depth of 2.0 m 
the central propeller induces, at the beginning of stern tunnels 
(x = -1.1m), a higher drop of pressure than side propellers 
(Fig.4). The results seem reasonable because the central pro- 
peller is closer to the waterway bottom than the side propel- 
lers. For smaller depths one may expect even a higher drop of 
bottom pressure and thus a greater destructive effect of run- 
ning propeller on waterway bottom. This regularity is illustra- 
ted in Fig.5 for the side propellers. The results of computations 
also showed that the rotational speed of propellers within the 
considered range had no effect on the distribution of bottom 
pressure. 


0.4 


Bottom pressure at h = 2.0 m , V, = 0.25 m/s 
0.2 


Cp [-] 


x [m] 
1 0 i 2 3 -4 -5 


—?— beneath side propellers, n = 30.0 rps 
—t— beneath central propeller, n = 12.5 rps 


Fig. 4. Distribution of bottom pressure at near zero ship speed (stern tran- 
som at x = 0; side propellers at x = - 0.318; central propeller at x = - 0.234) 


0.2 
Bottom pressure beneath side propellers 


V, = 0.125 m/s 


x [m] 


1 0 
—e h=20m, 


-1 -2 -3 -4 -5 
n = 30 rps —E— h=2.0m , n=20 rps 
—ġ— h=5.0m , n=30 rps 


Fig. 5. The effect of water depth and propeller 5 rotational speed on the 
distribution of bottom pressure beneath side propeller at near zero ship 
speed (stern transom at x = 0; side propellers at x = - 0.318) 


THE EFFECT OF CANAL CROSS-SECTION 
ON PROPELLER PERFORMANCE 


The computations of propeller performance in the canal of 
trapezoidal cross-section were carried out in the scale of 1:14. 
The slope of banks was 1: 4. The draught of barges was assu- 
med equal to 1.4 m. In the computations the barge train was 
positioned in various distances from the banks (Fig. 6). For the 
considered pushboat the computations were carried out sepa- 
rately for the central and side propellers as the interaction be- 
tween propellers was experimentally proven negligible. 


The results for ship scale are presented in Tab.2 and 3. The canal width bo is measured at the canal bottom. 
The advance coefficient J was calculated with account for the effective wake fraction We. 


90 m 


h = 5.0 m (0.357 m) 


4.5 m (0.321 m 


10 m(0.714 m 


25 m(1.786 m) 


Fig.6. Cross-section of the canal and considered positions of the barge train 


Tab. 2. Performance of the side propeller in the canal (the propeller 5 Tab. 3. Performance of the central propeller in the canal (propeller s 
diameter D = 0.690 m, and its rotational speed n = 828.5 rpm) diameter D = 1.100 m, and its rotational speed n = 345.2 rpm) 


Water depth h = 1.7 m Water depth h = 1.7 m 


Canal width bo 9m 


0.834 | 0.753 | 0.788 | 0.707 | 0.640 | 0.558 

fraction We i i i 
Thrust coeff. Kr| 0.331 | 0.319 | 0330 ! 0315 | 0.325 ! 0.303 

Efficiency | 0.050 ; 0.127 | 0.062 | 0.147 | 0.092 | 0.212 Efficiency n | 0.079 | 0.289 | 0.114 | 0310 | 0.155 | 0.362 
Advance coeff. J] 0.030 | 0.080 | 0.038 | 0.094 | 0.057 | 0.142 Advance coeff. J| 0.068 ' 0.289 | 0.081 | 0.247 | 0.112 | 0.302 


Water depth h = 2.0 m Water depth h = 2.0 m 


Canal width bo Canal width bo 20m 


Ship speed Vs ' Ship speed Vs 
Nominal wake i i i Nominal wake i i 
; 0.632 ! 0.588 | 0.672 ! 0.629 | 0.778 | 0.74 . 0.437 | 0.402 | 0.469 ' 0.419 | 0.564 ! 0.529 
fraction Wn t H A fraction Wn ` 1 
l 0.020 ' 0.234 | 0.043 | 0.245 | 0.117 ' 0.348 
fraction We i í 
Thrust coeff. Kr| 0.296 | 0.274 | 0.305 ' 0.293 | 0.311 | 0.304 Thrust coeff. Kr] 0.393 ' 0.330 | 0.416 | 0.365 | 0.421 + 0.381 
Efficiency | 0.198 | 0.283 | 0.184 ' 0.254 | 0.158 | 0.203 Efficiency 7 | 0.294 ' 0.422 | 0.278 : 0.407 | 0.261 ' 0.371 
Advance coeff. J| 0.127 | 0.188 | 0.122 ! 0.175 | 0.103 ! 0.135 Advance coeff. J} 0.232 ' 0.362 | 0.226 : 0.362 | 0.208 : 0.314 


Water depth h = 3.6 m Water depth h = 3.6 m 


Canal width bo 50 m Canal width bo 
Ship speed Vs i i i Ship speed Vs 
1.5 3.0 1.5 3.0 1.5 3.0 [m/s] 
0.222 | 0.216 | 0.250 | 0.244 | 0.302 ' 0.296 

fraction wn ! i i 
- 0.285 | - 0.004 | -0.279 : 0.009 | -0.23 ' 0.060 

fraction We ! ! ! 
Thrust coeff. Kr] 0.265 : 0.217 | 0.282 : 0.255 | 0.284 + 0.259 

Efficiency» | 0.293 : 0.417 | 0.279 + 0.400 | 0.270 | 0.383 


Advance coeff. J| 0.202 + 0.316 | 0.201 + 0.312 | 0.193 + 0.295 Advance coeff. J| 0.288 | 0.470 
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Selected results are presented in Fig. 7 through 12. 

The bottom pressure at the position where slope begins is 
presented in Fig. 7 and 8. At a lower water depth and lower 
canal width the destructive effect of moving ship over water- 
way bottom is stronger. 


Bottom pressure at h = 1.7 m , V, = 1.5 m/s 


0.5 


—e Canal 1,b) =9m —l— Canal 2 , bo = 20m 


—»-— Canal 3 , bọ = 50m 
Fig.7. Longitudinal distribution of bottom pressure at the position 
where slope begins, for various values of canal width 


Bottom pressure at bọ = 9 m, V, = 3 m/s 


+> h=17m = h=2m — h=36m h=5m 


Fig.8. Longitudinal distribution of bottom pressure at the position 
where slope begins, for various values of water depth 

The effect of canal width and water depth on performance 
of side propellers is presented in Fig. 9 through 12. The results 
show that the effect of banks is pronounced and depending on 
water depth. At the lowest depth h = 1.7 m the propeller thrust 
decreases when canal width increases. At the higher depths 
(h = 2.0 m and h = 3.6 m) the trend is opposite. An increase 
of water depth at a constant canal width causes a drop of 
propeller thrust. This effect can be explained by changes in 
water inflow to propeller that is reflected in variation of wake 
fraction. 


22990 
ty Go Ww GO 
-NUA 


0 10 20 30 40 50 60 
—eh=17m,V,=15 m/s —m—h=17m,V,=3 m/s 


—*—h=2.0m,V,=1.5 m/s —*—h = 2.0 m , V, =3 m/s 
Fig.9. The effect of canal width on thrust generated by side propeller 
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Thrust of side propeller at h = 3.6 m and h = 5.0 m 
0.29 


0.28 
0.27 
0.26 
_ 0.25 
~~ 0.24 
20.23 
0.22 
0.21 
0.20 
0 


bg Im] 


10 20 30 40 50 60 


—@ h=3.6m,V,=15m/s =- h=36m,V,=3 m/s 
—#— h=5.0m,V,=3 m/s 


Fig.10. The effect of canal width on thrust generated by side propeller 


Thrust of side propeller 


0.34 
0,32 
0.30 
= 0.28 
¥ 0.26 
0.24 
0.22 
0.20 


1 2 3 4 him] 5 
=% b= 9m, V,=1.5 m/s ==- b= 9m, V,=3 m/s 
—t— bọ=20m, V, =1.5 m/s = bo) =20m, V, =3 m/s 
Fig.11. The effect of water depth on thrust generated by side propeller 

Thrust of side propeller 


1 1.5 2 
—e b=50m, V, = 1.5 m/s 
—t— bo) =50m,V,=3.0 m/s 


3.5 4 45 5 
—& bo) =20m,V,=1.5 m/s 
=é bo) =20m,V,=3.0 m/s 


2.5 3 


Fig.12. The effect of water depth on thrust generated by side propeller 


In the case of the central propeller the effects of canal width 
and water depth are similar to those above presented. 


If the ship is not positioned in the centreline of the canal the 
side propellers operate in different conditions. At the same ro- 
tational speed the starboard and portside propellers generate 
different thrust and the ship tends to turn if the turning moment 
is not compensated by rudders. 


CONCLUSIONS 


O It was proved that the applied computational method is 
a valuable tool to predict the performance of ducted pro- 
pellers in shallow water, even at zero ship speed. 


O Both the results of the model tests and computations reveal 
that the effect of water depth on thrust and pull at zero ship 
speed is insignificant. 


O 


In the case of the considered barge train (dẹ = 1.4 m and 
dp = 0.6 m) the computed drop of bottom pressure beneath 
the running propeller is several times less than drop of pres- 
sure beneath bow and stern of barges. In general, it is the 
matter of ship and propeller arrangement as this factor le- 
ads to a greater destruction. 


The effect of canal banks and bottom on propeller perfor- 
mance in the case of the considered vessel is evident at the 
water depth below 3.6 m (h/d, = 2.6 and h/d, = 6.0) and the 
canal width below 20 m (bọ/B = 2.2). Propeller thrust at the 
lowest depth and width is 17% to 25% higher than that in 
wide and deep canal, depending on ship speed. Unfortuna- 
tely the propeller efficiency drops then dramatically. 


NOMENCLATURE 
bo - canal width 
B - moulded beam P-P 
Cp - pressure coefficient (= ———=*) 
d  - moulded draught —pVè 
D - propeller diameter 
h  - water depth 
J  - advance coefficient T 
Kr - thrust coefficient (=—;—7) 
Loa - overall length of ship P” D 
Lw, - length at waterline 
n  - rotational speed of propeller 
P  - local static pressure 
Pee - reference pressure (static pressure far upstream) 
PS - plane of symetry 
T  - thrust of propeller 
V, - ship speed 
We - effective wake fraction 
W, - nominal wake fraction 
x  - longitudinal co-ordinate, measured from propeller disc 
n  - propeller efficiency 
p  - water density 
V - volumetric displacement 
Indices 
b - of barge 
p - of pushboat 
Acronyms 
CFD - computational fluid dynamics 
INBAT - innovative barge trains for effective 
transport on inland shallow waters 
RTD - research and technological development 
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so C onference - 


A scientific seminar 
at Technical University 
of Koszalin 


On 24 March 2005 the first-in- this-year scientific 
seminar of the Regional Group of the Section on Explo- 
itation Foundations, Machine Building Committee, Po- 
lish Academy of Sciences, took place at the Department 
of Mechanical Engineering, Technical University of Ko- 
szalin. 

According to the announced program of the seminar it 
was devoted to role of heuristics in didactic and scientific 
research activity. 4 papers on the topic were presented by 
scientific workers of the University, namely : 


æ Heuristics in creative solving the educational problems 
by J. Markul 

% Heuristics in innovation and development of products — 
— exemplified by a didactic task — by J. Plichta 

# Heuristics as a method of solving logistic problems in 
manufacturing processes — by G. Jurkowski 

% Heuristics applied to solve converse tasks in machine 
building — by W. Tarnowski. 


After interesting discussion the seminar's participants 
visited laboratories of the Department. Additionally, the 
hosts informed about very broad spectrum of didactic and 
scientific research lines of the Department's activity cove- 
ring : mechanics and machine building, new technologies, 
solid body physics, materials engineering, optimization, 
automation, robotics and control, production process au- 
tomation, mechatronics, data processing, modelling, simu- 
lation, artificial intelligence methods, biotechnologies, 
biochemistry, food processing, design, bionics. 

An important form of the activity is the developing of 
cooperation with foreign scientific centres. In present the 
Department maintains contacts with three German and 

ae French centres. 
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High-cycle fatigue criterion for anisotropic 
metals under multiaxial constant 
and periodic loads 


Janusz Kolenda 
Gdańsk University of Technology 


ABSTRACT 


Periodic stress with Cartesian components given in the form of Fourier series is conside- 
red. To account for the mean stress effect the modified Soderberg s formula is employed. 
An equivalent stress with synchronous components is defined. The design criterion 
7 for a finite fatigue life of metal elements is formulated. It covers the conditions of both 
static strength and fatigue safety in the high-cycle regime and includes material constants 
which have simple physical interpretation, can be determined by uniaxial tests, are related 

directly to the applied loads, and can reflect material anisotropy. 


Key words : design criteria, multiaxial loading, periodic stress, mean stress effect 


INTRODUCTION 


In static problems, the load capacity of metal elements is 
referred to the yield strength and/or ultimate strength. In the 
case of multiaxial static loads, a reduced uniaxial stress, equi- 
valent to the original one in terms of effort of the material, is 
taken into account. For ductile metals such stress model can be 
defined, for instance, with the aid of the Huber-von Mises- 
-Hencky distortion-energy strength theory [1,2]. In the case of 
dynamic loading conditions, the design criteria of engineering 
elements made of ductile metals are frequently based on their 
fatigue limits and/or S-N curves (Wöhler curves) [3+5] and, 
if the stress state is multiaxial and the stress components are 
proportional to each other, also on the distortion-energy theory 
[6,7]. In [8] the design criterion for an infinite fatigue life of 
metal elements under the stress with non-proportional compo- 
nents has been formulated with the aid of the average-distor- 
tion-energy strength hypothesis [9] and theory of energy trans- 
formation systems [10]. In this paper the design criteria for 
a finite fatigue life of metal elements under the stress with non- 
-proportional components in the high-cycle regime are consi- 
dered. 


A CRITERION OF FINITE FATIGUE LIFE 
UNDER UNIAXIAL STRESS 


Ifa metal element is subjected 
to an axial load producing the stress : 


6, (t)=ctasinat (1) 
where : 


a -stress amplitude 
c -mean value 
œ -circular frequency, 
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in design for an infinite fatigue life various „failure diagrams” 
or equations can be used [5,11]. For example, in [8] the So- 
derberg’s equation has been utilized. Its necessary modifica- 
tion to indicate finite fatigue lives resulted in the following 


formula [11] : 
a=0 (l-c/R,) (2) 


where : 
Re - tensile yield strength 
o - amplitude of the fully reversed stress 
at a given number, N, of cycles to failure 
a - amplitude of the stress (1) which leads 
to that fatigue life. 


For the relation between o and N in the high-cycle fatigue 
regime the following equation of the S - N curve is commonly 
accepted : i 

No =K (3) 
where : 


K -fatigue strength coefficient 
m -fatigue strength exponent. 


Introducing the safety factor : 


f= N = T (4) 
N Ty 
where : 


Na = Tg /(27)- required number of stress cycles 
to achieve a given design life Tg 
T- time to failure under the stress (1) 


and partial safety factors : 


su feh (5) 
N, c 


Mee (6) 
a 
one gets 
f=f,(1-f,')™ (7) 
The criterion in question reads : 
f>1 (8) 
that is : 
f,d-f,)" 21 (9) 
Eq. (9) can be rewritten as : 
f +f <1 (10) 
1e 
1/m 
Ni åt us c<l (11) 
K R, 


It is seen that Eq.(11) covers both the condition of static 
strength and the fatigue safety requirement. This advantage will 
be maintained in what follows. To be on the safe side, in the 
case of c <0 it is recommended to insert into Eq. (11) the com- 
pressive yield strength Rec< 0 in place of R, . 


A CRITERION OF FINITE FATIGUE LIFE 
UNDER OUT-OF-PHASE STRESS 


On the basis of the average-distortion-energy strength 
hypothesis [9], the stress with Cartesian components : 


0; (t) =c;+a, sin(wt + B; ) (12) 
i =X, y, Z, XY, YZ, ZX 
can be modelled by the reduced stress [12] : 
O,,(t) =c,,+4,, sin wt (13) 


equivalent in terms of effort of the material under the stress 
with the components (12), where : 


ci, aj, Bi- mean value, amplitude and phase angle 

of i-th stress component, respectively 

- mean value and amplitude of the reduced 
stress, given by (for the sake of brevity the 


stress components Oz , Oyz and 67x have 


Ceq » deq 


been dropped) 
ew) 2 \1/2 
Ceg = (Cx +C,- CC, +3c,, ) (14) 
2.2 2 41/2 
aq =[a, +a}—a,a, cos(B,-B,)+3a,,]~ (15) 
Consequently, the criterion (11) becomes : 
1/m 
N 1 
ag ty, = 1 (16) 
K R, 


or, in explicit form 
N 1/m 
Tia? 2 41/2 
(e Jisi+al—aa, eos. B,) +305, + 


l 
toites ec, +304)? <1 (17) 


e 


In particular, the criterion (17) may be useful 
for isotropic steels where [1] : 


=——R 
es J3 e t 
where : 


Res - shear yield strength 
Fp - fatigue limit under fully reversed bending 
F; - fatigue limit under fully reversed torsion. 


With Eqs (18), the criterion of finite fatigue life for steel 
elements under out-of-phase stress can be also expressed as : 


N 1/m 
ca [ay+a,—a,a, cos(B,-B,)+ 


+E ja? J? + (19) 
F|? 
t 
1 R, ) 
+—[c} +c? -ce +| — Je J? <1 
R, Ras 


A CRITERION OF FINITE FATIGUE LIFE 
UNDER OUT- OF-PHASE STRESS 
FOR ANISOTROPIC METALS 


Metals may be isotropic or anisotropic with respect to any 
type of behaviour, such as thermal conductivity, electrical resis- 
tivity, thermal expansion, plastic deformation, or fatigue 
strength. This paper is concerned with differences in yield 
strengths and fatigue limits in different directions as well as 
under loads of different modes. For instance, uniform stress 
field under axial load versus nonuniform one under bending 
contributes to significant difference in the relevant fatigue li- 
mits. 

Anisotropy can be introduced by cold-working operations 
such as rolling, forging, drawing or stretching. For example, 
even small amounts of strain can cause a considerable diffe- 
rence in the longitudinal and transverse yield strengths. In 
other words, as a result of the longitudinal prestrain the tensi- 
le and compressive yield strengths are different in the trans- 
verse direction, as well as in longitudinal direction, though to 
a smaller degree [13]. Cold rolling of a wide plate and cold 
drawing of a tube over a mandrel both result in a lengthening 
and thinning, with no change in the transverse dimension 
(width or diameter), and so the resulting anisotropy is similar 
in the two products. In both cases, the anisotropy exhibits 
three mutually perpendicular planes of symmetry. In drawing 
of a rod (or tubing without a mandrel), the two transverse 
strains are equal, and so the anisotropy reveals an axis of sym- 
metry. 

For different metals and different fabrication methods the 
degree of anisotropy varies within wide limits; some metals 
are so nearly isotropic that it may be difficult to detect a diffe- 
rence in properties in different directions, such as in the case of 
the yield strengths in the normal and rolling directions of hot- 
rolled mild-steel plate; however other metals may be highly 
anisotropic. 

In order to take into account the material anisotropy and 
load modes, in this paper the following modification of Eq. 
(19) is suggested : 
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NV". T fe.) fe F rr YU 
— =a, |+| ay |-| =a, || =a, Jcos(B -B + =|| —a„ | |+ 
K F, F, FJ F F, || F, 
2 2 2 i (20) 
R, e R, R, R, R.,, 
+— c, |+ c, |- (N c, |+ Ča <1 
R,|| R, R, R, “|R, Ra || Ry 
where : 


F, - fatigue limit under fully reversed load relevant to the normal stress component of the amplitude a, 

Fxy- fatigue limit under fully reversed load associated with the shear stress component of the amplitude ayy 
R, - absolute value of the yield strength relevant to the constant load inducing the normal stress component c, 
R,y- yield strength associated with the constant load giving the shear stress component cy, . 


The remaining material constants Rj and F; are defined analogously. 


Eq. (20) yields the requested criterion : 
1/2 


ifi 2 Ve 2 

N a; C, cC 

ca F, > Aa F F ay cos(B.— P,) |+ > sa = <1 (21) 
x" Sy 


i i i i 


A CRITERION OF FINITE FATIGUE LIFE UNDER MULTIAXIAL CONSTANT 
AND PERIODIC STRESSES FOR ANISOTROPIC METALS 


Let us consider a periodic stress with Cartesian components given by Fourier series : 


0,(t)=c, + > a sin(po,t +B) 1=X, y, XY (22) 
p=l 
where : 
(A - mean value of i-th stress component 
Sips Bip - amplitude and phase angle of p-th term in 


Fourier expansion of i-th stress component 
@ọ = 27/To - fundamental circular frequency 
To - common period of the stress components. 


Following the results obtained in [8] and [12], the stress components (22) are modelled by the equivalent stress components : 


oD A= +a sin(w,,t + @;) 1=X, y, xy (23) 
where [8] : 
R 1/2 
c = c; ad = Fa aca a cos(@,— P, )= Dan m cos(B,,— B) (24) 
p=l 
and [12]: 
O = KO, k = Round(k) (25) 
1/2 
2 (pap) + oap) Z Èr a spå yp cos(B,; Po) i 391 Paa 
m pe (26) 


ae t Le Zana yp cos(B,.- Bp) +3193, 
p=! 


Here k is a natural number, qj, is the phase angle of i-th equivalent stress component and eq is the equivalent circular 
frequency. Eq. (26) has been derived by means of the theory of energy transformation systems [10] under assumption that the 
material is ductile and isotropic. Its modification for ductile and isotropic materials, similar to that of Eq. (17), leads to : 


z È -De poos BaBy) 


kK=| 7 i (27) 
5) > D son p 008(Bap- Pal 
i p=l F p=! FF 
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Adaptation of Eq. (21) to the equivalent stress with the components (23) gives : 


2 


1/2 2 


1/2 
1/m A (eq)„ (eq) (eq) (eq) a (eq) 
Na a (ev a a c; C, € 
—i | F i — — _¢os(9_- + i -> 7 <1 (28) 
K ° py F, FEF, we > R; R,R, 
where : 
y 2al (29) 
d 
2T 


Hence the criterion of finite fatigue life for engineering details made of anisotropic ductile materials 
and subjected to multiaxial constant and periodic loads becomes : 


T l/m 2 1/2 2 1/2 
O acf a 1 < c. Cx 
eq sd Ip i xy 
F — — ) a „a „ cos(B,— + - —— |<l (30) 
2nK : py 2, F] EE >, py Cos(B,— By) py Ri] RR, 
EXAMPLE With Eq. (21) one obtains : 
Assumptions 


Suppose the stress components 
and material data (in MPa) are : 


0, (t) =c,+a, sin@t — normal stress component 
O(t) =c,,+a,, sin @t — shear stress component 
c,=50 a,=90 c,=40 a,=60 
R, =260(steel25) R,, =160 
R. =310 — yield strength in bending 
F,=Z,,=180 F=Z,, =110 


F, =Z. =150 — fatigue limit under 
fully reversed tension - compression. 


Task 
Compare the criteria (17) and (21) in the following cases : 


A. the normal stress component is caused by bending moment 
B. the normal stress component is caused by axial force. 


Solution 


A common and conservative strategy in classical fatigue 
design for applying the distortion-energy strength theory is to 
separately compute the equivalent stress corresponding to the 
mean stress and the equivalent stress corresponding to the am- 
plitude [14]. Similar results have been obtained with the aid of 
the average - distortion - energy strength hypothesis, Eqs. (14) 
and (15), which will be used in the considered Example. 


A. Eqs (14) and (15) yield : 
Coq = (CZ +307, )"? = (50° +3-40°)"? = 85.44 MPa 


aq =(a, +3a,,)” =(90° +3-60°)"? =137.48MPa 
So, Eq. (17) becomes : 


1/m 
137.48 Mi + aa 
K 260 


Hence : 


1/m 
Z) < 48.83.107 
K 


N,\™ (90? 602) (502. 402) 
— |180 zt >| + zt z|<1l 
K 180° 110 310° 160 

Le., 


1/m 
a < §2,.74:10 * 
K 


B. In this case Eq. (17) leads to the same results 
as with the bending stress, that is : 


1/m 
G < 48.83.107 
K 


whereas Eq. (21) gives : 

1/m 2 2 N12 2 2 
Nu) 130/20. 4 0°), (50, 40 
K 150° 110 260° 160 


1/m 
ca < 46.90-10* 
K 


1/2 


The differences between these results are self-explanatory. 


CONCLUSIONS 


@ The finite fatigue life design criterion covering the condi- 


tions of both static strength and fatigue safety of metal ele- 
ments under multiaxial constant and periodic loads, has been 
formulated. 


@ The presented criterion includes material constants which 


= have simple physical interpretation 
= can be determined by uniaxial tests 
=> are directly related to the applied load 


= can reflect material anisotropy. 


Similarly as in vibration problems, it is not possible to state 


a definite upper limit to p in (22), (24), (26), (27) and (30), 
since this depends upon the nature and origin of the con- 
sidered stress variations, as well as upon the influence of 
minor terms on fatigue endurance of various materials, 
which may be the source of uncertainties. 
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NOMENCLATURE 


- stress amplitude 

- amplitude of i-th stress component (i = x, y, Z, XY, YZ, ZX) 

- amplitude of the reduced stress 

- amplitude of p-th term in Fourier expansion of i-th stress 
component 


ai - amplitude of i-th component of the equivalent stress 
c - mean stress value 

Cj - mean value of i-th stress component 

Ceq - mean value of the reduced stress 

oe - mean value of i-th component of the equivalent stress 
f - safety factor 

fa» fs - partial safety factors 


, F, - fatigue limits under fully reversed bending and torsion, 
respectively 


F; - fatigue limit under fully reversed load associated with the 
stress component of the amplitude a; 
k - natural number obtained by rounding the number K 
K - fatigue strength coefficient in equation of the S-N curve 
for tension-compression 
m - fatigue strength exponent in equation of the S-N curve 
for tension- compression 
N,N, - numbers of zero mean stress cycles to cause failure at the 
stress amplitudes o and a, respectively 
Na - required number of stress cycles to achieve a given 
design life 
Re, Res- tensile and shear yield strengths, respectively 
Ri - yield strength associated with the constant load inducing 
the stress component c; 
t - time 
T - time to fatigue failure 
Ta - design life 
To - stress period 
Bi - phase angle of i-th stress component 
Bip - phase angle of p-th term in Fourier expansion of i-th 
stress component 
K - quantity given by Eq. (27) 
(o) - stress amplitude satisfying Eq. (3) 
Oa - stress produced by axial load 
0; - i-th stress component 
Seq - reduced stress 
of - i-th component of the equivalent stress 
Q; - phase angle of i-th component of the equivalent stress 
0) - circular frequency 
Mo - fundamental circular frequency of the periodic stress 
Weq -~ equivalent circular frequency 
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a C onrerence SZ, 


Safety at sea and marine 
environment protection 


These were the topics of 9 Technical Scientific Con- 
ference held in Kołobrzeg on Baltic Sea coast on 2 + 3 
June 2005. The conference program contained 17 papers 
presented during two sessions : 


Safety of navigation and rescue of lives at sea 
(13 papers) 
Marine environment protection (4 papers). 


Out of the presented topics, 8 papers dealt with local 
problems, and general ones were given in the following 
papers : 


> Safety of maritime turism in the polar regions, with 
Antarctica as an example — by J. Frydecki and A. Wol- 
ski (Martime University of Szczecin) 

> Development trends in modelling the areas of search 
at sea — by Z. Burciu (Gdynia Maritime University) 

> Contemporary methods for classification of free- drift 
phenomenon — by M. Drogosiewicz and A. Wojcik 
(Polish Naval University) 

> Problems associated with dynamical determining the 
search area — by T. Budny (Gdynia Maritime Univer- 
sity) 

> Introduction of manoeuvrability standards in the low- 
-speed range as an element of improvement of port 
operation safety — by T. Abramowicz-Gerigk (Gdynia 
Maritime University) 

> Modelling the life- raft sinkage probability 
by L. Smolarek (Gdynia Maritime University) 

> Probability of radar detection of rigid inflatable boats 
(RIB) fitted with active radar reflectors — by A . Szklar- 
ski (Gdynia Maritime University) 

> Application of ultrafiltration technique to bilge water 
purification — by Z. Jozwiak and A. Kozłowski (Mari- 
time University of Szczecin) 

> Research on deoiling process of oil/water emulsion by 
means of ceramic bulkheads — by J. M. Gutteter — Gru- 

~~ dzinski (Maritime University of Szczecin). p 
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Forces exciting vibrations 
of ship’s hull and superstructure 


Lech Murawski 
Ship Design and Research Centre, Gdansk 


ABSTRACT 


The paper presents identification of main excitation forces of vibrations of ship hull and its 
superstructure, as well as investigation of influence of the analyzed excitation on ship 
vibration level. The following excitation forces were analyzed: dynamic pressure on the 
transom deck, propeller-generated hydrodynamic forces, as well as internal and external 
engine-generated forces. Rigidity and attenuation characteristics of lubricating oil film 
were taken into account in the analysis of the excitations transmitted by the propulsion 
system bearings (stern tube bearing, intermediate bearing, engine 5 main and thrust bea- 


rings). The dynamic characteristics of ship's power transmission system were also taken into account. 
Results of the analyzes were verified by comparing results of the computations and measurements. 


Keywords : hull vibration, superstructure vibration, main engine vibrations, excitation forces, 
pressure pulses, propulsion system’s reactions, crankshaft-propeller phasing 


INTRODUCTION 


The following main forces exciting vibrations of ship hull 
and superstructure (Fig.1) were considered [1,12] : 


1. propeller-induced pressure pulses upon the ship transom 
deck 


2. dynamic reactions of the thrust bearing to non-coupled lon- 
gitudinal vibration of the power transmission system (ge- 
nerated by axial hydrodynamic force of the propeller) 


3. dynamic reactions of the thrust bearing to coupled longitu- 
dinal vibration of the power transmission system (genera- 
ted by kinematic couplings of the crankshaft and hydrody- 
namic couplings of the propeller) 


4. dynamic reactions of the radial bearings (stern tube bearing, 
intermediate bearings and main bearings of the engine), 
generated by hydrodynamic radial forces and moments in- 
duced on the propeller 


5. non-balanced moments of the main engine, generated by 
gas and inertia forces of piston-crank system. 


The analysis was performed for the cargo ships (mainly for 
medium size containers) with typical propulsion system — a con- 
stant pitch propeller directly driven by a low speed main engine. 
In the paper an example analysis was made for a 2700 TEU 
container carrier. The ship particulars were as follows: ship 
length — 207 m, maximum deadweight — 35600 t, speed — 
— 22.7 knots. The ship was powered by MAN B&W 8S70MC-C 
engine of 24840 kW rated output at 91 rpm. The five-blade 


propeller was of 7.6 m diameter and 42400 kg in weight. All 
conclusions are valid for the ship loaded with containers from 
1200 TEU up to 4500 TEU in number. For the analysis of the 
ship’s hull and superstructure vibrations Patran-Nastran softwa- 
re was used, whereas the analysis of excitation forces (dyna- 
mic calculations of the power transmission system) was made 
mainly with the use of the author’s DSLW FEM - based soft- 
ware dealing with coupled torsional-bending-longitudinal vi- 
brations of power transmission system [7] and DRG one con- 
cerning shaft line whirling vibrations [8]. 


Fig. 1. Excitation forces of ships hull and superstructure vibrations 


POLISH MARITIME RESEARCH, No 4/2005 15 


MARINE ENGINEERING 


Reactions of stern bearing, shaft line bearing and main bea- 
rings of the engine are generated by bending vibration of the 
shaftline. The reaction level and distribution depend on shaft line 
alignment [5,11] and on stiffness-damping characteristics of lu- 
bricating oil film in the bearings [6,8]. The stiffness-damping 
characteristics of the oil film were found with the use of the au- 
thor’s BRG software based on Finite Difference Method [6]. 
Dynamic rigidities of ship hull structure within areas of bearing 
foundations were determined with the use of Nastran software. 

Analysis of bending vibration of the propeller shaftline sho- 
wed a low level of bending vibration, nevertheless the level of 
dynamic reactions of the bearings (particularly the stern tube 
bearing) was significant. It may be a source of excessive vibra- 
tions of ship hull and superstructure. 

Reactions of the thrust bearing are generated by longitudi- 
nal vibration of the power transmission system. [3,7]. Two ba- 
sic sources of longitudinal vibration excitation forces may be 
distinguished : excitations generated by the main engine [4,7] 
(coupled torsional-bending-longitudinal vibrations) and exci- 
tations generated by the propeller (coupled and non-coupled 
hydrodynamic forces). During the analysis total reactions of 
the thrust bearing, generated by all kinds of kinematic couplings 
were found. These reactions result from longitudinal vibrations 
of the power transmission system, coupled with torsional-ben- 
ding vibrations [3,7]. Non-linear characteristics of thrust bea- 
ring’s oil film were taken into account during the calculation. 

Reaction phases generated by non-coupled vibrations were 
determined with respect to propeller blade (phase ,,zero” was 
assumed for the upper position of the blade). Reaction phases 
generated by coupled vibrations were determined with respect 
to the first crank of the engine (phase ,,zero” was assumed for 
the top dead centre). The amount of total reactions (generated 
by coupled and non-coupled vibrations) depends on mutual 
phasing of propeller blades and crankshaft. With a suitable ar- 
rangement of angular position of the propeller (in relation to 
the crankshaft) the resultant vector of the reactions can be re- 
duced to algebraic difference of modules of the summed reac- 
tions. The optimum setting angle of the propeller with respect 
to main engine cranks depends on rotational speed of the pro- 
pulsion system, minimized harmonic component and load con- 
dition of the ship. 

The reaction of longitudinal vibration damper should be 
applied to the axis of the crankshaft due to its full symmetry. 
To find the application point of thrust bearing reaction is more 
difficult. The thrust bearing (of Mitchell type) consists of pads 
supported at their edges and arranged at the perimeter of thrust 
disk. Generally the pads do not fill completely the full perime- 
ter, being installed at the lower part of the thrust disk. In this 
case the axis of thrust bearing reaction should be lowered. 
Another solution consists in applying the thrust bearing reac- 
tion in line with the crankshaft axis and adding a suitable ben- 
ding momentas a transverse dislocation of force does not chan- 
ge the loading system if adequate moment is added. The effect 
of neglecting this moment on the vibration level of ship hull 
and its superstructure should be investigated during further 
numerical analysis. 

Regardless of the bearing forces, vibrations of ship hull and 
its superstructure are generated by pressure pulses induced on 
the plating above the propeller (transom), as well as by unba- 
lanced internal forces of the main engine. The pressures, their 
distribution and phase shift are specified by designers of the 
propeller. The unbalanced forces and moments of the main 
engine are the last significant forces to be analyzed. In contem- 
porary main engines all the external forces are generally ba- 
lanced. The remaining moments which excite the following 
vibrations of engine body, should be investigated: vertical-lon- 
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gitudinal (ML), vertical-horizontal (MH) and torsional vibra- 
tions around the vertical axis (MX). For the seh in question, 
(ie. MAN B&W 8S70MC-C), only the 1° 4 5™ and 8 
harmonic components of excitations are as from the 
point of view of ship hull and superstructure vibrations. 


Two main structural models of the 2700 TEU 
container carrier were used for the computations : 


¢ the model for the ship in the design load condition (Fig.2) 
and that for the ship in the ballast condition [9,10]. 


Each of the models 
contained more than 42000 degrees of freedom. 


Fig. 2. Structural model of 2700 TEU container carrier 
in the design load condition 


Results of the computations are presented for the represen- 
tative group of the points (all the points, except the bridge wings, 
are placed in the ship plane of symmetry) located as follows : 


BOW- at hull bow on the main deck 

DTR -at deck transom edge on the main deck 
SBF -at fore bottom of the superstructure 
SBB —at aft bottom of the superstructure 
MEF -at main engine fore cylinder heads 
MEB - at main engine aft cylinder heads 

STL -at bridge wing (left) 

STF -at fore top of the superstructure 

STB -at aft top of the superstructure 


KINDS OF EXCITATION 


LULL 


The analyses of forced vibrations were made for the contai- 
ner carrier in the design load condition (Fig. 1) without added 
mass of water. The wet model of ship hull was used to find the 
total forced vibrations (generated by complex sources of exci- 
tations) and to compare them with measurement results. Com- 
putations for the container carrier in the ballast condition were 
also carried out because it was the condition of the ship during 
the measurements. 

The propeller induces a variable field of water pressure exci- 
ting vibration of ship transom deck. The vibrations are transmit- 
ted to all parts of the ship structure. Vibrations excited by the 
first-blade (5" harmonic) component of propeller-induced exci- 
tations were analyzed. The image of forced vibration velocities 
of the ship hull, superstructure and main engine is shown in Fig.3. 

Hydrodynamic forces and moments are induced by rota- 
ting propeller. Only the variable hydrodynamic forces acting 
in line with propeller shaft and generating non-coupled axial 
vibration of the propulsion system are discussed. Phases of these 
vibrations are strictly associated with the propeller. The non- 
-coupled axial vibrations of the power transmission system ge- 
nerate excitations to ship hull and superstructure through the 
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Fig. 3. Vibration velocities of the container carrier, 
excited by water pressure field on the transom deck 


thrust bearing and axial vibration damper. The forced vibra- 
tion of ship hull, superstructure and main engine is presented 
in Fig.4. 


Fig. 4. Vibration velocities of the container carrier, 
forced by axial vibrations of its propulsion system 


Torsional-bending vibrations are generated by ship propul- 
sion system. Bending of crankshaft cranks and torsional vibra- 
tions of the crankshaft and propeller generate coupled longitu- 
dinal vibrations of the power transmission system. The longitu- 
dinal vibrations induce dynamic reactions on the thrust bearing 
and axial vibration damper, generating excitations to ship hull 
and superstructure. Phases of the coupled longitudinal vibra- 
tions depend mainly on the position of main engine crankshaft. 

The coupled longitudinal vibrations are mainly due to gas 
and inertia forces generated in the engine cylinders. Therefore 
arelatively wide spectrum of harmonic components may be of 
significant importance. In the considered case the first eight 
harmonic reactions (8-cylinder engine) are significant. There- 
fore eight variants of computation of the forced vibrations of 
ship hull and superstructure were carried out, i.e. for all signi- 
ficant harmonic components of exciting forces. Distribution of 
forced vibration velocities for the ship is similar to that shown 
in Fig.4. 


The general image of forced vibration amplitudes 
(velocities) for all the considered harmonic components 
is shown in Fig.5. 
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Fig. 5. Longitudinal vibration velocities 
forced by coupled longitudinal vibrations 


As in the case of other kind of exciting forces, the highest 
vibration level appears in the upper part of ship superstructure. 
The 5" harmonic component is dominant as in the case of exci- 
tations generated by the propeller. The levels of expected vi- 
brations generated by the coupled and non-coupled longitudi- 
nal vibrations of the propulsion system are similar. Phases of 
vibrations generated by the non-coupled vibrations depend on 
a position of the propeller blades, whereas phases of vibrations 
generated by the coupled vibrations depend on a position of 
the crankshaft cranks. Therefore a mutual angular position of 
the propeller and crankshaft may significantly affect vibration 
level of ship superstructure. This problem is further investiga- 
ted below. Side harmonic components (with respect to 5™ com- 

` th th ae A 
ponent), i.e. 6 and4 components may be dominating for vi- 
brations of the main engine body and ship hull. 

The transverse forces and moments induced by rotating 
propeller generate bending vibration of the shaftline. The ben- 
ding vibrations generate dynamic reactions of radial bearings 
of the propulsion system, i.e. the stern tube bearing, interme- 
diate bearing and main bearings of the main engine. The resul- 
ting bending vibrations of the shaftline generate forced vibra- 
tion of ship hull and superstructure. 

Fig.6 presents the vibration velocities of ship hull, super- 
structure and main engine, forced by transverse vibration of 
the propulsion system. 


Fig. 6. Vibration velocities of the container carrier, 
forced by transverse vibration of the propulsion system 


POLISH MARITIME RESEARCH, No 4/2005 17 


MARINE ENGINEERING 


MARINE ENGINEERING 


Non-balanced external dynamic forces and moments appear 
as aresult of operation principle of piston engines. In ship en- 
gines all external forces and high-order moments are balanced. 
However some of the moments remain still non-balanced. The 
moments generate longitudinal, transverse and torsional vibra- 
tions of the engine, transmitted to other structural parts of the 
ship. Values of the non-balanced moments are given in the en- 
gine documentation. For the considered engine the following 
harmonic components are significant: 1, 3, 4, 5, and 8. 

Velocities of forced vibrations of the ship for: 1harmonic 
(excitation of L and X type), 3™ harmonic (excitation of X type), 
5" harmonic (excitation of X type) and 8" harmonic compo- 
nent (excitation of H type) for rated rotational speed of the 
propulsion system are shown in Figs.7 to 10. 

During computation of the ship vibrations forced by the 
external non-balanced moments of the main engine, the exci- 
tation moments of the order 1., 3., 5. and 8. appear significant 
for the considered type of the main engine (8S70MC-C). The 
1% harmonic component excites mainly vibration of the hull 


Fig. 7. Vibrations forced by 1%” harmonic component (L and X type) 
of main engine 5 external moments 


Fig. 8. Vibrations forced by 3"! harmonic component (X type) 
of main engine ï external moments 


Fig. 9. Vibrations forced by 5" harmonic component (X type) 


of main engine ï external moments 
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Fig. 10. Vibrations forced by git harmonic component (H type) 
of main engine s external moments 


structure, whereas the g” component is responsible for vibra- 
tion level of the superstructure and main engine (particularly 
in transverse direction). The excitation forces generated by non- 
-balanced moments of the main engine are responsible parti- 
cularly for increasing the level of transverse vibrations. For the 
remaining orientations of vibration other excitation forces are 
dominating. 


OPTIMIZATION OF MUTUAL ANGULAR 
SETTING OF PROPELLER 
AND CRANKSHAFT 


Two sources of excitations should be taken into account 
when analyzing the dynamic behaviour of ship hull and super- 
structure: 


+ the excitations generated by the propeller and 
+ those generated by the main engine. 


For the propeller-induced excitations the 5™ harmonic com- 
ponent is dominating. The engine-induced excitations are repre- 
sented by full spectrum of harmonic components (8 harmonic 
component is main for the engine in question), but also the 5! 
harmonic component is here significant. According to the abo- 
ve presented analyses the levels of expected vibrations forced 
by the both excitation sources are similar. Phases of the vibra- 
tions generated by the propeller depend on a position of the 
propeller blades, whereas the phases of the vibrations genera- 
ted by the main engine depend on a position of crankshaft cranks. 
Therefore mutual phasing of the propeller and crankshaft may 
considerably affect vibration level of ship superstructure. 

Ship superstructure vibrations ofeach orientation (longitu- 
dinal - x, transverse - y, vertical - z) may be optimized (to redu- 
ce its level to a minimum). Fig.11 shows the optimum relative 
angles of angular position of the propeller blade against the 
first crank of the crankshaft. 

Only one angular position of the propeller can be chosen. 
The highest vibration level is expected for longitudinal orien- 
tation of vibration (in x direction). The assumed optimum an- 
gle œ = -20.3° gives the minimum level of longitudinal vibra- 
tion at rated rotational speed of the engine. 

Fig. 12 presents the total vibration level at the assumed 
mutual angular phasing of the propeller and crankshaft. It also 
presents the maximum vibration level (for the worst angular 
phasing of the propeller : & =+15.7°), as well as the minimum 
one. 


Optimal angle 


Phase [deg] 


60 65 70 75 80 85 90 95 
Main engine revolution [rpm] 


Fig. 11. Optimum phasing angles of the propeller 
against the main engine crankshaft 


x-direction optimization 


Velocity [mm/s] 


60 65 70 75 80 85 90 95 
Main engine revolution [rpm] 


Fig. 12. Longitudinal vibration velocities of the superstructure (bridge 
wing) for the optimum phasing of the propeller and main engine crankshaft 


Minimum level of longitudinal vibration may be obtained 
within a wide range of engine rotational speeds for the opti- 
mum angle of propeller angular phasing with respect to the 
crankshaft. It could be achieved because the optimum angles 
for this orientation of vibration (x direction) were almost con- 
stant in function of engine rotational speed (Fig.11). At the 
same time the vibration levels for other orientations of vibra- 
tion (e.g. transverse) may be not optimal, for instance vibra- 
tion velocity values may be higher than the minimum ones 
theoretically possible. 

The optimum angle value (œ = -20.3°) is near to that obta- 
ined from another method based on analysis of vibration pha- 
ses of deckhouse (œ = -23.8°). This angle was found by opti- 
mizing phases of the forces exciting coupled and non-coupled 
longitudinal vibrations of the power transmission system. Chan- 
ges in phases of structural response and other excitation sour- 
ces were not taken into account. The correctness of the opti- 
mum angle was verified by the measurements carried out on 
the ship in question for different angular positions of the pro- 
peller. During sea trial an unacceptable vibration level of the 


ship superstructure was measured. After changing the relative 
propeller angle the vibration level measured during next sea 
trial was reduced as much as three times. 


SUMMARY VIBRATIONS — 
— VERIFICATION BY MEASUREMENTS 


The above presented analyzes are based on the hull structu- 
ral model without added mass of water, for the ship in design 
load conditions. In order to verify the presented method the 
forced vibrations should be computed both for the ballast and 
load conditions by using the mathematical model of the ship 
with added mass of water. In such mathematical model all the 
formerly analyzed kinds of excitations were assumed to act 
simultaneously. Additional independent computations were 
carried out for the vibrations forced by the main engine and 
then for those forced by the propeller. The results of these com- 
putations can be directly compared with the measurement re- 
sults [2]. 


The velocities of the forced vibrations of the ship at rated 
rotational speed of the propulsion system are shown in Fig.13. 


Fig. 13. Total vibration velocities of the container carrier 
in design load condition 


Figs.14 and 15 present the diagrams of longitudinal vibra- 
tion velocities of the superstructure for both considered ship’s 
load conditions, separately for different types of excitation. 
Possible minimum and maximum vibration levels are also 
shown. The following notation was used: Prop —excitation by 
propeller, M.E. — excitation by main engine, Max —maximum 
vibrations, Min — minimum vibrations, SUM — expected sum- 
mary vibrations. 


ò Design load (x-direction) 
1 T T T T 


Velocity [mm/s] 


60 65 70 75 80 85 90 95 
Main engine revolution [rpm] 


Fig. 14. Total velocities of longitudinal vibrations of bridge wing 
(ship in design load conditions) 
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Ballast load (x-direction) 


Velocity [mm/s] 


60 65 70 75 80 85 90 95 
Main engine revolution [rpm] 


Fig. 15. Total velocities of longitudinal vibrations of bridge wing 
(ship in ballast conditions) 


Figs. 16 and 17 present a comparison of the computation 
results and ship trial measurement results. The results were 
converted into rms values of vibration velocity to obtain the 
customary presentation. The comparison is presented for two 
essential reference areas : the bridge wing of the superstructure 
and the transom deck. 


STL (x-direction) 


V - rms [mm/s] 


60 65 70 75 80 85 90 95 
Main engine revolution [rpm] 
Fig. 16. Verification of the bridge wing vibration: computation results ver- 
sus measurement results 

The forced vibration level is clearly higher for the ship in 
ballast conditions when compared with that for design load 
conditions. Sea trials are usually carried out in ballast condi- 
tions of the ship and therefore even a correctly designed ship 
may show excessive vibrations of the superstructure during 
service. Vertical - horizontal vibrations of the upper part of 
superstructure are dominating. Comparable vertical vibrations 
may be also observed for the transom deck. 

Magnitude of excitation forces generated by the main engi- 
ne is similar to that generated by the propeller. Therefore, when 
analyzing the hull and superstructure vibrations, at least the 
excitation forces generated by pressure pulses on the transom 
deck and the forces generated by longitudinal vibration of the 
power transmission system should be taken into account. Ships 
of an identical hull structure may show a high or relatively low 
level of vibration depending on whether its propulsion system 
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is suitably designed. Therefore the complete dynamic analysis 
of the ship should be carried out simultaneously with its de- 
sign process — such analysis should be one of its elements. 


DTR (z-direction) 


V - rms [mm/s] 


60 65 70 75 80 85 90 95 
Main engine revolution [rpm] 
Fig. 17. Verification of the transom deck vibration: computation results 
versus measurement results 

The verification presented in Figs 16 and 17 showing a sa- 
tisfactory compliance of the numerical computations with me- 
asurement results, confirmed that the applied assumptions and 
computation methods were correct. 


RECAPITULATION 


O The excitation forces generated by non-balanced moments 
of the main engine are responsible particularly for increa- 
sing the level of transverse vibrations. In the case of dyna- 
mic analysis of the container carrier in question other exci- 
tation forces are dominating (in longitudinal-vertical direc- 
tion). 


O If levels of the total vibration generated by the propeller 
and main engine are similar to each other the vibration le- 
vels may be optimized by changing the relative angular 
position of propeller blades and crankshaft. Vibrations of 
each orientation may be optimized. The longitudinal vibra- 
tions expected to be ofa high level (excited by main excita- 
tion sources: pressure pulses on deck transom and thrust 
bearing reaction) are usually optimized. 


O Forced vibration level is clearly higher for the considered 
ship in ballast conditions than in design load conditions. 
The level of excitations generated by the main engine is 
similar to those generated by the propeller. Therefore, when 
analyzing the hull and superstructure vibrations, at least the 
excitation forces generated by pressure pulses on the tran- 
som deck and the forces generated by axial vibration of the 
power transmission system should be taken into account. 


O Ships of an identical structure may show a high or relative- 
ly low level of vibration depending on whether their pro- 
pulsion system is suitably designed or not. Therefore a com- 
plete dynamic analysis of the ship should be carried out in 
line with its design process, being one of its elements. 


O Correctness of the used assumptions and computation me- 
thods was confirmed by the measurements (Figs 16 and 
17) as the presented results of the numerical computations 
and measurements appeared to be in a satisfactory com- 
pliance. 


ACRONYMS 


BRG -bearing 

DSLW - coupled vibration of shaftline 
DRG - bending vibrations 

FEM - Finite Element Method 
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Numerical simulation of heat flow processes 


Andrzej Mielewczyk 
Gdynia Maritime University 


ABSTRACT 


Numerical simulation method is based on mathematical models of physical processes, which 
describe run of a given process with a various accuracy. This paper presents a method for 
elaborating the computer simulation models based on differential equations. This makes it 
possible to obtain a high accuracy of representing dynamic features of a real process. The 
derived results are given in the form of numerical series. An example of modelling the 
dynamics of heat transfer through a flat wall is presented in the second part of the paper. 
Basing on the model one can simulate the process of operation of ship engine as well as 
such auxiliary devices as a cooler, evaporator, condenser, tank heating system etc. 


Keywords : simulation, dynamic model, heat transfer, Z - transform 


INTRODUCTION 


Computer simulation has presently become a common tool. 
Simulators are used in didactic process, scientific research, and 
training courses, especially for operators. This way they sub- 
stitute real technical devices, for economical and safety reasons. 
Today operational and maintenance conditions of machines are 
closer and closer to those represented in simulators. Educatio- 
nal process of to-be operators of machines starts from using 
the simulators. 

The simulation makes it possible to represent difficult and dan- 
gerous operational situations and to learn to give an appropriate 
response to hazards. A boundary between a real situation and 
simulated one becomes obliterated. A simulation must be close 
to reality in order an operator always could consider a given 
situation as real one. By continuous improving the simulators 
a higher and higher level of education and safety can be ensured. 

Developments in process automation and control have made 
it necessary to have at one’s disposal a model of a given pro- 
cess [8, 9, 12]. 

By comparing a model of a given process with its real run 
it is possible to determine an instantaneous point of operation 
select appropriate control parameters and generate a correct 
message for operator, that rises safety level of the process. Si- 
mulation has also found wide application in training of sea- 
-going ship’s operators, both navigators and marine engineers. 
Hence simulators have become an indispensable instrumenta- 
tion of maritime academies. 

Simulation approach has also its important place in for- 
ming new scientific theories and widening the knowledge. In 
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this way one looks for confirmation of laws in engineering, 
physics, economy and medicine. An expected chance for spa- 
ce missions is also checked by means of simulations. 

In this paper has been presented the problem of modelling 
heat flow processes by using Z - transform applied to diffe- 
rential equation of a modeled process. This facilitates elabo- 
rating the digital models on which the work of any simulator 
is based. 


SIMULATION MODELLING 


During process simulation one tends to exactly represent 
statical and dynamical features of a given process. A simula- 
tion model should appropriately respond to external signals, 
and in the case of some kinds of simulators it should provide 
real-time responses. Real processes are complex and depen- 
ding on many parameters. Only the simulation models based 
on differential equations are capable in providing a high con- 
formance with real processes [10,14]. To find solutions of such 
equations is time-consuming. Modelling the objects of complex 
parameters leads to partial differential equations [1, 3, 4, 11]. 
For this reason one often introduces simplifications in expense 
of accuracy of process representation. For example, either sim- 
ple algebraic equations are introduced, or process statical and 
dynamical features are separately considered . However the so 
seperated dynamical term does not account for all variables of 
the object in question. 

Applying another approach one substitutes finite difference 
equations for differential equations, that impairs accuracy of 
solution. In solving partial differential equations by means of 


the finite difference method a condition for stability of solution 
is introduced [13]. For instance the heat transfer equation : 


dU(x,t) _ r 0°U(x,t) 
ot Ox? 


obtains, after digitization respective to time and spatial varia- 
ble, the following form of response in successive time instants : 


(1) 


UG, tas )= = aa Ux, ? ta )+ 


T T 
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To reach a higher calculation accuracy a digitizing step, 
both in the time and space domain, is made shorter and simul- 
taneously the following stability condition has to be satisfied : 


1 

a <— 
Ax? 2 6) 
The time step T is imposed by the assumed digitizing step 
Ax. The application of a variable digitizing step is necessary 
to accelerate a response being always of recurrent type. How 
to find a compromise between simulation accuracy and its du- 
ration time is an open question. During simulation process 
much time is usually devoted to graphical presentation of 

a model. 


This author has proposed to combine the advantages 
of solving differential equations and finite difference ones 
in order to determine process models. 


The digitization is applied respective to the time variable t 
only. Equations in such a form can be solved with the use of 
the transform Z. The time variable t passes into z and becomes 
a parameter. The so obtained equation can be solved respective 
to spatial variables as continuous one by integrating it and ap- 
plying the known methods for solving partial differential equa- 
tions, e.g. the method of separation of variables or that of suc- 
cessive integral transformations [5]. In this phase dynamic re- 
lations of modeled process is not solved. Now the inverse trans- 
formation Z of the initial function is made to obtain its course 
respective to time. 

The proposed approach makes it possible to obtain a solu- 
tion equivalent to continuous one. This way the digitization 
in the domain of spatial variables as well as the necessity to 
satisfy the associated stability condition, is avoided. The di- 
gitization in the time domain only enables to choose indpen- 
dently a time step limited only by an assumed solution accu- 
racy. The inverse transform Z is calculated by expanding the 
function into numerical series by means of the following for- 
mula [6]: 


ee I dp" (4) 

If an expansion of a given function into Taylor series re- 
spective to the variable p exists then it is always possible to 
determine a solution. This is much easier to do than to use the 
Laplace transform [2, 7, 10]. The final solution contains all 
process variables and the time variable in a discrete form. The 
result is yielded in the form of numerical series, which leads 
to simple procedures of numerical programming; the calcula- 
tion error depends on digitizing the function respective to time, 
only. 


In each computational step the process parameters take con- 
stant values. In Eq.(1) the parameter a is constant and not de- 
pendent on time and position. The parameters can be changed 
in successive computational steps. The equation is reduced to 
a differential equation of constant coefficients, and a linear form 
of a given model is obtained. 


EXAMPLE 
HEAT TRANSFER THROUGH A WALL 


Heat transfer through a flat wall occurs in many ship sys- 
tems and technical devices. The phenomenon is described by 
Eq.(1) which can be solved by using the proposed method. 
After application of the transform Z the Eq. (1) takes the follo- 
wing form: 


WU(x,z 1 z-l 
(x )_ 1 zl 2-0 (5) 
dx aT z 
where : 
À 
a = — 
pc 
After transformation one can obtain a parabolic 2™4 order 


differential equation of constant coefficients. Its solution is as 
follows : 
(6) 


U(x,z)=Ae™+Be™ 


where : 


r, >= 
i27 7 


In order to solve Eq.(5) the following 
boundary conditions have been defined : 


(7) 


1. Temperature on the inner side of the wall varies 
in compliance with the assumed function : 


U(0,z)=A+B=U,(z) (8) 


2. The outer side of the wall is insulated, hence : 
dU(x = w,z 
= dU(x = w,z) _ )= 0 (9) 
dx 


The solution of Eq.(5) is described 
by the following formula : 


—M(Are™ + Bre’ 


enw) 4 an(x-w) 


U(x,z) =U, (z) = (10) 


e7" +e 


The wall’s outer side temperature 
is described by the formula : 


) 2 
w {z-l w jz-l 
e vaT E, e vaT = 
Finally, it is necessary to calculate 
the inverse transform Z by using the relation (4) : 


Ge. glee 
=e '* +e '% 


where : 


U(w, z) =U,(z (11) 


(12) 


WwW 


vaT 
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The inverse transform is obtained this way : 


n-l n-2 
on 3) 22 1 (n-k-1)(2n-2k-3)! okt! 
f(n =0,-1, Do 3 a 2 2 knl(2k—1)\(n—2k—1)) sinh & + 
n-3 n-2 (13) 
o 2° 2 = _ 
+ 2cosho,0 272 (2n 2k ~3)! at cosh 


rm = KS nl(2k +1)!(n - 2k -2)! 
A choice of k value is realized in accordance with the following principle : 
k= 1 for odd values ofn ; k=2 for even values ofn 


The exciting function Ug(z) may be arbitrary, however this leads to the operation of convolution. If Up(z) is assumed a step 
function then the response is a sum of terms of given series. This way the mentioned operation of convolution can be avoided. 
For realization of other functions the approximation by a staircase function is advised. Such approximation method is used for 
digital control algorithms. 

In the time domain Eq.(11) takes the form : 
U(w,nT)=U, zi (14) 
ay ta; +a, ta; +... taata, 


where : 
_ o ; 1. 1 , 3 : 3.4 l ,. 
ag = 2cosh a ; aj =-Qsinha ; a, =-—asinhat+—a cosha ; a, =-—asinh a +— a" cosh æ -— a” sinh a 
4 4 24 24 24 


n-l n-2 
a, =- 1- Pua) T PL - e 1){2n = 2k — 3} aa sinha + 
27°? ni(n—2)! ao kn! (2k -1)!(n-2k-1)! 
n-3 n-2 
+ yy nas 3 CES a: o.**? cosha (1a) 
= 2 n!(2k +1)!(n-2k—2)! 
After dividing by the series a, the response at discrete instants is obtained. For the exciting step signal it amounts to : 
Es 3 n= 3) a tanha 
2°" nl(n-2)! 
asa 7 
U(w,nT)=U, : + y 2 5 — a n-k-I)n-2k-3) gap, , tanh?" a+ (16) 
cosha| £ 442" kn!(2k-1)!(n-2k-1)! i 
n-3 n-2 
: k 
pa __On=2k-3)) _gauap tanh a 
=e n!(2k +1)!(n — 2k — 2 
The coefficients Dm are derived from the recurrence algorithm 
based on the Pascal triangle principle, namely from the formulae (17 and 18) : 
k=1 1 
k=2 1 -2! 
k=3 1 -3! 
k=4 1 -10-2! 4! 
k=5 1 -10-3! 5! 
k=6 1 -91-2! 35-4! -6! 
k=7 1 -91-3! 35-5! -7! 
k=8 1 -820-2! 966-4! -84-6! 8! 
k=9 1 - 820-3! 966-5! - 84-7! 9! 
k=10 1 - 7381-2! 24970-4! - 5082-6! 165-8! -10! 
k=11 1 - 7381-3! 24970-5! - 5082-7! 165-9! -11! 
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D(2k +1,m)=(-1)" ID, (2k -1,m)+ (2m +1) D,(2k-1 m)|(2m +1)! p 
D(2k,m)= C1)" D, (2k -2,m)+ (2m +1} D, (2k -2,m)(2m)! 
k i! NOMENCLATURE 
m . . 
Dym = > (- 1) k! (i z k)! Dkk, i) (18) a — coefficient of temperature equalization [m?/s] 
i=m “Al : a,  — coefficients of numerical series 
For m=0 the Euler numbers are yielded [15] : A, B — integration constants 
k c  — specific heat of a material [J/kgk] 
= Aa D — coefficients of expanded form of a initial function 
Dyo > D{k,i) E (19) n  — natural number 
: i20 . Tı, fy — zero loci of a partial differential equation 
The successive values of the function are as follows : t  — times] 
— T -— digitizing period [s] 
U(w,0) 0 1 U _— temperature function 
= U, — initial values 
U(w,T ) ~ U, Tee w  — wall thickness [m] 
cosh & 
i 1 x — spatial variable 
Z, p — arguments of a function in the complex variable domain 
U(w,2T)= U, h l+ 2 a tanh & à  — thermal conductance [W/mK] 

cosh Q p  — density [kg/m3] 
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Below in the figure results of the example calculations 
performed by using the described method are presented. 
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Dynamic characteristics of heat transfer through the wall 
of 0.1 m in thickness, the excitation U(t) = 1[deg], for the following 
materials : 1 — copper; 2 — aluminium; 3 — brass; 4 — steel. 
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ding the preceding ones. The so obtained solution is equi- 
valent to continuous one. 
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An algorithm for determining permissible 
control inputs to unmanned Underwater Robotic 
Vehicle (URV) fitted with azimuth propellers 


Jerzy Garus 
Polish Naval University 


ABSTRACT 


The paper deals with synthesis of automatic control system for an unmanned underwater 
robotic vehicle. The problem of determining permissible propulsive forces and moments 
necessary for optimum power distribution within a propulsion system composed of azi- 
muth propellers (rotative ones). To allocate thrusts the unconstrained optimization me- 
thod making it possible to obtain a minimum-norm solution, was applied. A method was 
presented for assessing propulsion system capability to generate propulsive forces (set 
control inputs). For the case of lack of such capability an algorithm was proposed making 


modification of their values and determination of feasible propulsive forces (i.e. permissible control 
inputs), possible. A numerical example which confirmed correctness and effectiveness of the proposed 
approach, was also attached. 


Keywords : underwater vehicle, propulsion system, azimuth propeller, control 


INTRODUCTION 


Underwater Robotic Vehicles (URV) play an important role 
among various technical means used for searching seas and 
oceans. The unmanned floating units fitted with propulsors and 
capable of maneouvring are designed to realize tasks at the 
water depth from a dozen or so to several thousand meters. 


Robot’s motion of six degrees of freedom is described 
by means of the following vectors [1, 3] : 


n=k, y, z, >, 6 w]’ 
v=lu, v, w, p, q, r|" (1) 
t=[X, Y, Z, K, M, NĪ" 
where : 


n- vector of location and orientation in an inertial re- 
ference system 

v — vector of linear and angular velocities in a hull- 
-fixed reference system 

t — vector of forces and moments acting on the robot 
in a hull-fixed reference system. 


Contemporary underwater robots are often and often equip- 
ped with the automatic control systems which make it possible 
to execute complex maneuvers and operations without any in- 
tervention of operator. The main modules of such control sys- 
tem are shown in Fig. 1. Its crucial element is the autopilot which 
—on the basis of comparison of a current location of controlled 
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object with its set values — determines forces and moments to 
be generated by the propulsion system so as to make the ob- 
ject’s behavior complying with that assumed. The thrust vec- 
tor corresponding with them is computed in the thrust distri- 
bution module and sent to the propulsion system as a control 
quantity. 


Disturbances 


Power 
distribution 
module 


Propulsion 
system 


Autopilot 


Fig. 1. Schematic diagram of a control system for underwater robot 


The control laws implemented in the autopilot , which make 
it possible to determine propulsive forces and moments, have 
a general character; they do not take into account constraints 
put on maximum and minimum values of the thrusts which 
can be developed by particular propellers. It can cause that the 
obtained solution would be unfeasible for the propulsion sys- 
tem. Such situation may worsen control quality and cause that 
robot’s behavior would greatly differ from that assumed. 


Therefore a two-phase procedure 
of power distribution is proposed (Fig.2). 


In the first phase propulsion system capability of genera- 
ting the set control inputs Tq would be assessed and the permis- 
sible control inputs Tj (i.e. such values of forces and mo- 
ments which the propulsion system is able to generate) would 
be determined. Their values would be so calculated as — ensu- 
ring operation of propellers to a saturation limit at most — not 
to cause a drastic perturbation in the robot’s motion control 
process. 


Block 
of allocation 
of thrusts 
to propellers 


Block 
of determination 


of permissible 
control inputs 


Fig. 2. Schematic diagram of power distribution module 


In the second phase on the basis of T} the thrust vector f 


would be calculated, i.e. the allocation of thrusts to particular 
propellers would be performed. 


THRUST ALLOCATION PROCEDURE 
FOR HORIZONTAL MOTION 


The solution used in majority of conventional unmanned 
underwater robots is a structure having longitudinal and trans- 
verse metacentric stability which ensures motion with small 
trim and heel angles. Hence the basic motion of such objects is 
their translation in horizontal plane at changing depth of im- 
mersion, being a motion of four degrees of freedom. 


It makes it possible to split the propulsion system 
into two independent subsystems, namely : 


> vertical motion subsystem (in vertical plane) 
> horizontal motion subsystem (in horizontal plane). 


The first produces the propulsive force acting along verti- 
cal axis, and the second ensures translational motion along lon- 
gitudinal and transverse axes and rotation around vertical axis. 

In this paper an underwater robot is considered fitted with 
the propulsion system having the arrangement of propellers 
shown in Fig.3. Its vertical motion subsystem consists of two 
vertically arranged, ducted screw propellers. To assess capabi- 
lity of the subsystem to generate the set force Zq is not a diffi- 
cult task as its absolute value cannot exceed the sum of maxi- 
mum thrusts developed by the propellers [4]. 


Fig. 3. Arrangement of the propulsion system fitted with six propellers 


The horizontal motion subsystem consists of four spatially 
arranged azimuth propellers producing forces along longitudi- 
nal and transverse axes as well as moment of force in vertical 


axis. To assess capability of the subsystem to generate the set 
forces Xq and Yq as well as the moment Ng is a complex task 
as each of the propellers contributes to generating both propul- 
sive forces and propulsive moment. Therefore it is necessary 
to have a procedure making it possible to assess whether the 
set control inputs are feasible, and in the case if the subsystem 
is not capable to produce them to modify them in such a way 
as to ensure their feasible values. Hence further considerations 
are limited to plane horizontal motion of the robot. 


Let Tq = [Ta , Ta2 , taz]! = [X4 ; Ya, Nal! stand for the 
vector of set propulsive forces and propulsive moment, 
and f= [f] , fo , f3 , f4]! - for the vector of thrusts developed by 
the propellers. 


Moreover let the components of the vectors 
be constrained by the following constraints : 


t - (cm)? <0 for j=13 (2) 
f - (es |? < 0 for i=1,4 (3) 


resulting from design parameters, arrangement and orientation 
of the propellers within the hull of the robot. 


The vector of forces and moment, T , is associated with the 
thrust vector f by means of the following relationship [1,2] : 


Ta = Tia) (4) 
where : 
T(Q) — arrangement matrix of propellers : 
T(a)= (5) 
cos(a, ) cos(a, ) cos(a, ) 
=| sin(a,) sin(a, ) sin(c,) 


d,sin(a,-@,) d,sin(a,-9,) d,sin(a, —@,) 

a= [0] , 0, ..., 4] — vector of thrust angles 

œ; — thrust angle of i-th propeller, i.e. the angle between ve- 
hicle’s longitudinal axis and direction of its thrust force 

Ọ; — orientation angle of i-th propeller, i.e. the angle between 
vehicle’s longitudinal axis and the line connecting the 
vehicle’s mass centre and the propeller axis 

d; — distance of i-th propeller from the vehicle’s mass centre. 


The quantities 0; , @j and d; are shown in Fig.4. 


Fig. 4. Arrangement of 4 propellers in horizontal motion subsystem 


The propellers are aimed at developing such thrusts which 
can ensure generating the vector of set control inputs tg = [T41 , 
Ta2 , Ta3]!. Allocation of thrust values to particular propellers 
is realized in the power distribution module. The problem of 
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determination of values the thrust vector f on the basis of the 
vector of set control inputs Tg is usually considered as an 
unconstrained optimization problem in which a minimum- 
norm solution is searched for. The solution has the following 
form [1,5]: f 

f=T (a)tg (6) 


where : 
Moore-Penrose’s pseudo-inverse matrix : 


T (a) =T'(o)[T@)- Toy" 


Its practical application is possible then and only then, when 
no demand of developing a thrust value exceeding the limit 
value (3) by anyone of the thrust propellers is declared. If it is 
the case then the set control inputs cannot be generated and 
a modification of their values is necessary, i.e. determination of 
the vector of permissible control inputs T} = [Th Tiz -Tis 
Together with the vector of control inputs the vector of azi- 
muth angles is computed. A way of determining their values is 
presented below. 


ALGORITHM FOR DETERMINING 
THE VECTOR OF PERMISSIBLE 
CONTROL INPUTS 


Let the robot’s propulsion subsystem ensuring its planar 
horizontal motion be consisted of four azimuth propellers of 
the following features : 


“* of the same type, hence the thrust f}"™ =f™ for 
i,k =1,4 
** located at the same distance from the mass centre, symme- 


trically against the robot’s plane of symmetry, namely : 


_ fmax 
= f 


T T == 
ta dandy ed S o =o for i,k =1,4 
% whose thrust angles satisfy the relationships : 
m at every instant ¢ for i,k =1,4 andi#k: 
T T 
a; (t) mod — = &,(t)mod — = a(t) 
2 2 
TU 
0< Ol min s a(t)< Ob max E 2 Ol nin 
(7) 


infa; (9| = sinjo, (t)]] = sinox(t)] 
lcos[a; (t)]| = |cos[o., (t)]] = cos[ox(t)] 


= |sin[o, (t) - 9, ]] = sin[a(t) + 9] 


= thrust angle change by the value + Aq occurs in all the 
propellers simultaneously. 


By virtue of the above given assumptions and the relation- 
ship (4) the forces and moment at the instant ¢ are described by 
the following equations : 


= $ coslo t 


(8) 
= cosa D sign {cos[a,(t)| tf, 
al- $sin fosh 
i (9) 
= sin[o wl sign {sin [a, (t)| W, 
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)-ọ; lf, = 


)-o, ]}f, 


a(t)]= Ky d; sin [a 
)+ old sign {sin lox (t 


where : 


(10) 
= dsin [a(t 


T[Q(t)] , T2[(t)]- propulsive forces along longitudinal 
and transverse axis, respectively 
73[0(t)]- force moment around vertical axis. 


(To make the mathematical description more clear, the 
time symbol ¢ in the notation of the thrust angle c(t) is fur- 
ther omitted). 


The maximum values of the quantities, possible 
to be generated by the propulsion system are as follows : 


Timax = max M (a) = 4 COS(O in Jf nax 
Qmin r S Op, 
Vmax eae i. (a) = 4 Sin (Ona lo el 1) 
T3 max J max 6 T; (x) T 4df nax 
for: 
ont 
7 Ọ 


An analysis of Eq. (10) shows that the propulsive moment 
T3(&) depends on the thrust angle œ through the term sin(& + @), 
where the angle @ is a time - independent design parameter. 


p 
Oo 


sin (a + p) 


x 
3 2 


T 
6 
Fig. 5. Influence of selected values of the angle ọ on sin(a + @) value 


In Fig.5 runs of the relationship sin(a + @) in function of 
the angle we (0, 2) are illustrated for some values of the angle 
@ contained in the interval (0,5-). They demonstrate that for 
a fixed value of the angle @ it is possible to determine such 
value of the propulsive moment which can be always genera- 
ted independently of a current value of the thrust angle a. The 
value is further marked T3max and calculated as follows : 


Tama < min t (a)= 
Adsin(Opint+ Pfa for QE 0,7) (12) 
~ ™ 1 
Ad si ff Ua 


Therefore, by applying the following constraint 
onto the propulsive moment Tq3 : 


Tis- Toa < 0 (13) 


its feasibility is ensured within the entire range 
of variability of the thrust angle a. 


Values of the propulsive forces which should be developed 
by the propellers to generate the moment T3max independently 
of a value of the thrust angle  , are determined by the relation- 
ship: 


- Tomax for QE o5) 
_]4d sin(c,,,,+ O) 4 (14) 
T3 max l T3 max for Ọ 2 x T 
4dsin(,,,. + 0) 4 2 


The reserving of a part of propeller’s output to generate the 
force moment T3max makes that the maximum values of the 
forces Timax and T2max Which can be generated by the propul- 
sion system, will be smaller. Their new values marked Tmax 
and Tomax, respectively, are now : 


Timax = 4 COS(O nin (Era fryma ) 
Tmax = 4 sin (Omax )(F ~ Fes ) 


max 
In further considerations it was assumed that : 


T 
Ta = [Tar ; Ta2 » Ta3] 
is the vector of set control inputs whose 
components satisfy the constraints : 


2 2 2 2 2 2 
Tar — Timax £ 0 2 Taz — Tmax £ 0 > Tas — Tmax = 0 


(15) 


and P is a point of planar coordinates (T42, Tq1). 


By analyzing the relationships (8) and (9) it was stated that 
the demanded propulsive forces Tg; and Tg2 can be simultane- 
ously generated by the propulsion system then and only then if 
the point P is located inside or at the edge of a geometrical 
figure shown in Fig.6. The figure is an asteroid described as 


follows : 2 2 2 
T ET =T 0 (16) 


whose vertices are the points of the coordinates : 


(T2max, 0) ; (0,Timax) ; (- Tomax, 9) 
and (0 , - Timax) respectively. 


4 


2 2 2 
Fig. 6. The asteroid described by the equation 7}, + Ti, = Trax 
and the position vector OP 


However if the point P lies outside the asteroid then the 
propulsion system is uncapable to develop set forces. Then the 
vector of permissible forces and moment t',=[t', , T! sys 
and the corresponding thrust angle «' should be determined. 


The schematic diagram of the algorithm for determining 
the permissible control inputs t} and the angle &' is shown in 
Fig.7. Input data for the algorithm are : 


> the set vector of control inputs Ty 
> the maximum value of the propulsive force Tj max 
> a current value of the thrust angle o 


and, the task of determining the permissible values of propul- 
sive forces T} and T!,, as well as of the thrust angle œ' is reali- 
zed for the following conditions : 


+ the propulsive moment is kept unchanged : T),=T,, 


7 T , f $ 
+ the mutual ratio of the forces ; ~+ = ~ is maintained. 


Ti Ta 


INPUT DATA 


Ta = [ta > Taz > Taal 
Ce 


Timax > 


CALCULATE: 


È , 1 
Tal > Tag > Q 


Ba d , 
Ta = [tay > Taz > Taal 


Fig. 7. Schematic diagram of the algorithm for determining permissible 
control inputs and thrust angle 
The algorithm in question consists of the following steps : 
1. Calculate the expression : 
2, 2 2 
3 3 r3 17 
Tå +T, — Tha LO (17) 


to check whether the point P = (Ta2 , Ta1) 
lies inside or at the edge of the asteroid. 


2. If the inequality (17) is true then apply substitutions : 


T a = Ta 
' = 

Taz = Ta, 
o'= 0 


and go to Step 4. 
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3. If the inequality (17) is false then calculate : 
3 


T Ta VP 2 
ta dl dl 
T dl sign(Ty, ) Le 4 + 1 T max 
a2 || \ Taz 
2 ME 
= 2 
oe ays dl 
Ta sign(t,, ) +1) Timax (18) 
d2 
a 
a'= =~ arctan} |! 


2 Taz 


4. Apply substitutions : 
TiS [th „Taz fal 
5. The end of the algorithm. 
The proof of Eq. (18) can be found in [6]. 


Numerical example 


Numerical calculations were carried out 
for the following data : 


ta =[700N ; -120N ; 50Nm]! 
T3min=50Nm , f nax =250N 
@=30° , d=0.4m 


The value of Timax was calculated for the worst case, i.e. 
for the angle Omin = 0° : 


Timax 7 4( fie fr) = 
=4 E ET. =750N 
4-0.4-sin(30°) 


Step 1 
Check if the point (Ta2 , Tq1) lies inside or at the edge 
of the asteroid, i.e. check the condition (17) : 
2 2 2 
Tir + Tay Tine SO 


2/3 


700° +(- 120)? — 7507 


= 20.6>0 
Step 2 


As the inequality has appeared false the point (Tg2 , Tq1) 
lies outside the asteroid. Calculate — by using (18) — permissi- 
ble values of forces and thrust angle : 


A. Calculation of the value of the force Th: 


2 


2 
f TiN 
T a2 > sign(T,, ) (3) + 1 Timax = 


d2 


2 E 
= sign(-120) (<a) + 750 =-85.9 
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B. Calculation of the value of the force T}, : 


2 
Tai : 
+ 1 Timax = 
Taz 


3 


Ta = sign(t,,) oa 
a2 


2 
=sign(700) 1° 750 = 501 


2 
700 \3 
120 


ed 
-120 


C. Calculation of the value of the angle a : 


3 
Tar 


Taz 


T 
a'=—~— arctan 
2 


Í 
008 gS irad=291" 
20 


T 
=— -arctan 
2 


Step 3 
Calculation of the vector of permissible control inputs : 
t =[501 ; -85.9 ; 50] 


To assess correctness of the calculations it was checked if 
the ratio of longitudinal and transverse forces was kept un- 
changed: 

Ta _ 700 on 


=—=-5.3 ; is 
Ta -120 Th, 


501 
-85.9 


= -5.83 


The obtained result showed that after correction the ratio 
of the forces was of the same value. 


RECAPITULATION 


O The presented paper concerns synthesis of an automatic 
control system for unmanned underwater robot, particular- 
ly its phase connected with power distribution in multi-pro- 
peller propulsion system. Planar horizontal motion of the 
robot executed by means of four azimuth propellers was 
considered. 


O For thrust allocation to the propellers an unconstrained 
optimization method was applied. To use the method prac- 
tically in a propulsion system of a limited power output 
a procedure which makes it possible to assess its capabi- 
lity of generating demanded forces and moment was 
elaborated. 


O For the case when the task appears unfeasible an algorithm 
was proposed allowing to modify propulsive forces by pro- 
portional decreasing their values and determining permis- 
sible ones, i.e. the forces which the propulsion system is 
capable to develop. 


O The performed numerical tests confirmed correctness of the 
applied approach. 


NOMENCLATURE 

d; — distance from i-th propeller to robot’s mass centre 

fi — i-th propeller’s thrust 

f — vector of thrusts 

K -— propulsive moment around longitudinal axis of hull-fixed 


reference system 


M -— propulsive moment around transverse axis of hull-fixed 
reference system 


N — propulsive moment around vertical axis of hull-fixed 
reference system 

p — angular velocity around longitudinal axis of hull-fixed 
reference system 

q — angular velocity around transverse axis of hull-fixed 
reference system 

r— angular velocity around vertical axis of hull-fixed reference 
system 

t  — time 

T(a)— arrangement matrix of propellers 

u — linear velocity along longitudinal axis of hull-fixed 
reference system 

v  — vector of linear and angular velocities in hull-fixed reference 
system 

w — linear velocity along vertical axis of hull-fixed reference 
system 

x — x- coordinate of vehicle’s position in inertial reference 
system 

X  — propulsive force along longitudinal axis in hull-fixed 
reference system 

y —y-coordinate of vehicle’s position in inertial reference 
system 

Y — propulsive force along transverse axis in hull-fixed reference 
system 

z — Z- coordinate of vehicle’s position in inertial reference 
system 

Z — propulsive force along vertical axis in hull-fixed reference 
system 


a; — i-th propeller thrust angle 

n — vector of location and orientation of vehicle in inertial 
reference system 

6 — heel angle 


v — linear velocity along transverse axis of hull-fixed reference 
system 

t  — vector of propulsive moments and forces 

Ta — vector of demanded propulsive moments and forces 

@; — i-th propeller orientation angle 

b — trim angle 

yw — course angle 
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jo — C onference -~ 
MECHANIKA 2005 


On 4 February 2005, 
Mechanical Faculty, 
Gdańsk University of Technology, 
organized the Scientific Conference : 


Mechanika 
on the occasion of 60 anniversary of its activity. 


Such conferences were earlier arranged in the years 
1995, 1997 and 1999. Their idea was to promote achieve- 
ments of scientific workers from mechanical faculties of 
the technical universities located in North Poland, i.e. in 
Gdansk, Gdynia, Szczecin, Koszalin, Bydgoszcz, Olsz- 
tyn, Białystok and Elbląg., as well as of the Institute of 
Fluid Flow Machinery, Polish Academy of Sciences, 
Gdańsk. It was mainly expected to draw interest of indu- 
strial enterprises of that region to those achievements. 
However effects of the attempts appeared unsatisfactory 
as technological parks, industrial fairs and branch confe- 
rences have become more effective occasions for contacts 
between scientific and industrial circles. 


As a result of the situation, this-year Conference was 
devoted to mutual presentation of selected research results 
and discussion on prospects of development of mechani- 
cal sciences within the frame of developing international 
cooperation. 


In compliance with the Conference’s program 35 papers 
were presented during three topical sessions : 


* Techniques and engineering processes 
of manufacturing (12 papers) 
* Drives and energy systems (12 papers) 
* Computer methods in mechanics (11 papers). 


Representatives of the Conference’s organizer - 
- with 17 papers - contributed the most 
to the elaboration of the presented papers. 


The remaining ones were prepared by authors from : 


the Faculty of Ocean Engineering and Ship Technology — 
— Gdansk University of Technology, Olsztyn Technical 
Agricultural Academy, Biatystok University of Techno- 
logy, Technical University of Szczecin, Warmia-Mazury 
University, Institute of Fluid-Flow Machinery of PAS — 
— Gdansk, Koszalin University of Technology, Polish Na- 
val University, the State Higher School of Engineering in 
Elblag, and ALSTOM company. 
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UNDERWATER TECHNOLOGY 


60 Anniversary 
Jubilee Academic Year 2004/2005 


The years 1945+2005 have been a period of dynamical 
development of Gdansk University of Technology which has 
today become the greatest technical university in the nor- 
thern region of Poland. The University may be proud of gra- 
duating over 80 thousand engineers in many disciplines, and 
of many achievements in the area of scientific research. The 
University’s employees have taken part in the activity of 
many international scientific associations and they have or- 
ganized many times international conferences giving impact 
to development in many domains of contemporary science. 
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Gdańsk University of Technology - main building (photo: C. Spigarski) 


These achievements justify undertaking a very compre- 
hensive, almost a forthnight lasting programme of celebra- 
tionof60 anniversary of that so much mertorious universi- 
ty. Beginning from the end of May this year many historical 
and scientific conferences, cultural events and light program- 
mes were held. It was also an exceptional occasion for arran- 
ging a gaudy-day of the university’s graduates to give them 
an opportunity to meeting in close circles of colleagues. 


Among nine today existing faculties of the University 
six of them have celebrated their own jubilee of 60 anni- 
versary because in 1945 they formed an origin for the fur- 
ther development of the University. 


The Faculty of Ocean Engineering and Ship Technology 
(formerly Faculty of Shipbuilding) also belonged to that gro- 
up. The Faculty much contributed to dynamical develop- 
ment of shipbuilding industry in Poland, which built more 
than 5000 ships of different types, among which a dozen or 
so was distinguished by international maritime bodies. 


Within the frame of the celebration programme of the 
Faculty’s jubilee, was held a solemn open session of the Fa- 
culty Council devoted to history of activity of the Faculty. 


> Education of many engineers and scientific workers with 
scientific degrees was acknowledged as the greatest achie- 
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vement of the Faculty. As many as 2226 B.Sc.diplomas 
and 3409 M.Sc. diplomas together were issued in the 
specialties of : Building of Sea-going Ships, Building of 
Inland Waterways Ships, Ship Power Plants and Machi- 
nes, Ship Equipment, Ship Technical Operation, Mana- 
gement and Marketing in Maritime Economy. Also, 329 
doctoral studies were completed. 


The 60 - year history of the Faculty has been described 
in the carefully edited book of 410 pages, titled : 
„Shipbuilding studies at Gdansk University 
of Technology in the years 1945=2005”. 


The Faculty’s graduates have been a basic, highly appre- 
ciated group of specialists employed in Polish industrial 
enterprises and shipping companies; many of them have 
carried out a fruitful didactic and scientific activity at the 
universities. 


A number of them have worked for foreign firms outsi- 
de Poland. Some of them, namely Mr. T. Berzowski, 
E. Haciski, T. Jaroszyński, W. Kossakowski, L. K. Ku- 
pras, R. Scheinberg and Z. K. Wójcik have imparted 
a piece of remembrance just on their proffessional acti- 
vity in various parts of the world. 


Finally, a solemn accent of the celebration was the cere- 
mony of giving the Faculty’s modern laboratory for fuel 
oils and lubricants the name of Prof. Janusz Staliński, 
deceased in 1985 former Rector of the University, Dean 
of the Faculty and multi-year Head of the Department of 
Ship Power Plants. Prof. Staliński gave a great impact to 
the Department’s development in respect to its didactic, 
scientific and design activities. He tightly cooperated both 
with industrial enterprises and shipping companies, that 
favourably influenced the course of studies in the field 
of ship power plants. 


Faculty of Ocean Engineering and Ship Technology 


